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We report first-principles calculations of local heating in nanoscale junctions formed by a single 
molecule and a gold point contact. Due to the lower current density and larger heat dissipation, the 
single molecule heats up less than the gold point contact. We also find that, at zero temperature, a 
threshold bias Vonset of about 6 mV and 11 mV for the molecule and the point contact, respectively, 
is required to excite the smallest vibrational mode and generate heat. The latter estimate is in 
very good agreement with recent experimental results on the same system. At a given external bias 
V below V onse t, heating becomes noticeable when the background temperature is on the order of 
~ e(V nset — V)/ks- Above Vonset, local heating increases dramatically with increasing bias but 
is also considerably suppressed by thermal dissipation into the electrodes. The results provide a 
microscopic picture of current-induced heat generation in atomic-scale structures. 
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Local heating occurs when electrons diffusing in a con- 
ductor release energy to the ions via scattering with 
phonons. The amount of heat generated in a given por- 
tion of the conductor depends on several factors: the 
strength of the clcctron-phonon interaction, the current 
density, the background temperature, and the electron 
mean free path. In a nanoscale junction, the electron 
mean free path is large compared to the dimensions of 
the junction. As a consequence, each electron releases 
only a small fraction of its energy during the time it 
spends in the junction. 0, 0, Q However, there can still 
be substantial local heating due to the large current 
density in nanoscale junctions. Indeed, such an effect 
has recently been observed in several atomic-scale struc- 
tures. |J,|5|,|£| Apart from its fundamental importance in 
solid state physics, the problem has thus gained renewed 
interest due to its possible impact on nanoscale electron- 
ics. At this length scale, all atoms of the junction and 
corresponding vibrational modes need to be treated ex- 
plicitely, and the electronic distribution calculated self- 
consistently with the correct scattering boundary condi- 
tions. Only a few theoretical investigations, mainly using 
the tight-binding approximation, have tackled this prob- 
lem at the atomic level. L3- LI 

In this letter, we derive a general expression for the 
power transferred by the electrons to the vibrational 
modes of a nanoscale junction and calculate the local 
heating using first-principles approaches. Combined with 
a model of heat transfer to the electrodes, the approach 
provides a microscopic and detailed picture of current- 
induced heat generation in nanoscale conductors. As an 
example, we discuss local heating for a single molecule 
between gold electrodes and a gold point contact. We 
find that the single-molecule junction heats up less than 
the gold point contact due to a larger heat dissipation 
into the electrodes and larger resistance to electron flow. 
In the case of the gold point contact, experimental data 
is available which is in very good agreement with the pre- 
dicted threshold bias, V onse t, for heat generation and cor- 
responding local temperature. At a given external bias 
V below V onse t, we find that local heating becomes no- 
ticeable when the background temperature is larger than 



~ eiVonset — V)/ks, but is also considerably suppressed 
by thermal dissipation into the electrodes. 

Let us start by considering a nanoscale structure be- 
tween two bulk electrodes. The many-body Hamiltonian 
of the system is 



where H e i 



H = H e i + H ion + H e i_ ion , 



is the electronic Hamiltonian, Hi, 
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5Xj Pf/2M i +X; i<j Vi „(Rj-Rj) is the ionic Hamilto- 
nian, and i? e /- l0 n = Vei-ioni^i - RjO describes the 
electron-ion interaction. Rj, Pj, and Mi are the coor- 
dinates, momentum, and mass, respectively, of the i-th 
ion. i"i is the coordinate of the i-th electron. In the adi- 
abatic approximation, the Hamiltonian Q has effective 

single-particle eigenvalues (r, Kii) with energy E 

and momentum Ky parallel to the electrode surface, cor- 
responding to electrons incident from the left (right) elec- 
trodes. We calculate the single-particle wavefunctions 
using a scattering approach within the density functional 
theory of many-electron systems. 

In order to calculate the electron-phonon interaction, 
we consider a small deviation Q; = R; — R? of the i- 
th ion from its equilibrium position R?. We introduce 
normal coordinates {<7i M } such that the /i-th component 
(ji = x, y, z) of Q. 4 is 
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The transformation matrix, A = satisfies 
the orthonormality relations: MiA^ju Ai^y v i — 

Sjv,j'v'- The Hamiltonian describing the ionic vibrations 
can then be decoupled into a set of independent harmonic 
oscillators: 
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i,fj,£vib. 



i,fj,£vib. 



where {^i^} are the normal mode frequencies and the 
summations are carried out for all normal modes. We 
calculate these modes from first-principles. 
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constant and has the form 
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where we have chosen to describe the electron-ion inter- 
action with the pseudopotentials V s (r. R,) for each i-th 
ion. hj v are the annihilation operators for the jV-th 
vibrational modes and satisfy the commutation relation 



The statistics of these modes are 



described by the Bose-Einstein distribution: 
(n jv ) = 1/ [exp (hw jv /k B T w ) - 1] 



(8) 



FIG. 1: Feynman diagrams of the four main electron-phonon 
scattering mechanisms contributing to local heating of the 
junction, (a) Cooling process due to absorption of a phonon 
from a left-moving electron, (b) Heating process due to the 
emission of a phonon from a left-moving electron, (c) and (d) 
are the equivalent mechanisms corresponding to the right- 
moving electrons. 



where T w is the local temperature of the junction. 

For each normal mode, we can now calculate the rates 
of thermal energy generated by the electron- vibration in- 
teractions corresponding to the first-order processes de- 
scribed in Fig.Q] These processes correspond to electrons 
incident from the right or left electrode that absorb (cool- 
ing) or emit (heating) energy because of the electron- 
vibration scattering. We evaluate these rates with the 
Fermi golden rule 



We now introduce the field operator, ^> = ^> L + 
which describes electrons incident from the left (L) and 
right (R) electrodes, where [j| 
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L(R)^L(R) 
E E 



r,K 
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Wff = 2nh(S k ,2 + (n jv )) I dE 



E 

iff 



^J" J E±Hu jv ,E 



— fE±nu ju )^ > E±tm : j U ^Ei 
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The coefficients a^ R ^ are the annihilation operators for 
electrons incident from the left (right) electrode. They 
satisfy the usual anti-commutation relations {a z E , a^L } = 
SijS (E — E'). We also assume that the electrons rapidly 
thermalize into the bulk electrodes so that their statis- 
tics are given by the equilibrium Fermi-Dirac distribu- 
tion, f^ R) = l/(exp[(£ - E FL{R) )/k B T e ] + 1), in the 
left (right) electrodes, where T e is the background tem- 
perature. 

We can now express (Qj)u m terms of boson annihi- 
lation and creation operators. Together with the field 
operator >F, the second-quantized form of H V ib and the 
electron- vibration interaction, H e i^ V ih, can be written as 
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where "R^± L" means interchange of labels R and L; the 
positive (negative) sign in Eq. J^J is for k — 1(2), corre- 
sponding to relaxation (excitation) of vibrational modes; 
is the partial density of states corresponding to 

L( R) 

fypf , whose sum is the total density of states. The 
term dk,2 corresponds to spontaneous emission. Finally, 
a factor of 2 due to spin degeneracy appears in Eq. 0. 

Since electrons can excite all possible energy levels of 
a mode with frequency Uj V , the statistical average (rij V ) 
is required. The total thermal power generated in the 
junction is therefore the sum over all vibrational modes 
for the four processes of Fig. JJ): 



P = 



E 

jis£vib. 
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and 



Hd-vib -EE E 

a,j3 E\ : E? ifi,ji/£vib. 



(6) 



where a, j3 — L,R ; i,j = 1, ...,JV; and jj,, v — 1,2,3. 
The quantity is the electron- vibration coupling 



Equation (|11|) is the central result of this paper. It al- 
lows for a first-principles calculation of the local temper- 
ature in a nanoscale junction. We now discuss results for 
two specific cases: local heating in a benzene-dithiolate 
molecular junction and a gold point contact (see schemat- 
ics in Fig. [2J . We also assume the left electrode to be 
positively biased. At zero temperature (i.e., T e = and 

T w = 0), the heating processes {w-jt' 2 }' corresponding 
to electrons incident from the left electrode, vanish due 
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FIG. 2: Local equilibrium temperature as a function of 
bias for a benzene-dithiolate molecular junction (dashed line, 
schematic in the right-lower corner) and a single-gold-atom 
point contact (solid line, schematic in the left-upper corner) 
at T e = K. No heat dissipation into the electrodes is taken 
into account. 



to the Pauli exclusion principle; the cooling processes 
{Wjv' 1 j an d {^jv' 1 }: corresponding to the transitions 



from a high to a low energy level of the modes, arc also 
prohibited because all modes are at the ground state. 
The only nonzero contribution to local heating is there- 
fore from W^' 2 , i.e. 
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A bias greater than V onse t — mhi{hujj U }/e is there- 
fore necessary to generate heat. The smallest vibrational 
frequency calculated from first principles || is about 6 
mV and 11 mV for the molecule and the point contact, 
respectively. The latter estimate is in very good agree- 
ment with recent experimental observations. When 
the bias is smaller than the threshold voltage, local heat- 
ing is only possible at nonzero background temperature 
T e . This is caused by a small fraction of thermally excited 
electrons which can induce level transitions of the nor- 
mal modes. The heating generated is substantial when 
k B T e « min{/iw JI/ } - (E FR - E FL ). 

At T e = and T w > 0, the junction heats up pro- 
gressively, and eventually an equilibrium temperature is 
reached when the heating processes (W^J 2 and W^J 2 ) 

balance the cooling processes (W^ 1 and W^J 1 ). The 
corresponding local equilibrium temperature is plotted 
in Fig. |21 for both the single molecule junction and the 



FIG. 3: The two modes contributing to heating at T e = and 
V < 10 mV. (a) The vibrational mode with energy hw ~ 6 
meV. (b) The mode with energy hw w 9 meV. 



gold point contact. It is evident from Fig.[21that the equi- 
librium temperature increases abruptly above the thresh- 
old voltage and it is already substantial at biases of only 
few mV. The reason for this dramatic increase resides 
again in the Pauli exclusion principle which suppresses 
the cooling processes considerably. 

In the case of the gold point contact, the longitudi- 
nal vibrational mode of the single gold atom (hw w 11.5 
meV) contributes the most to the local heating above 
V = 11.5 mV, even though the two transverse modes 
have slightly smaller frequencies (hw w 10.8 meV) and 
are excited at a lower bias. In linear response (V < 10 
mV), two modes contribute to local heating in the case 
of the single- molecule junction (see Fig. |3J). One, with 
energy hw 6 meV, corresponds to the central benzene 
ring moving up and down with respect to its equilibrium 
position (Fig. Ela)). The second, with energy hw « 9 
meV, corresponds to the central ring moving sideways 
with respect to its equilibrium position (Fig.^Jb)). The 
first mode contributes to local heating almost three or- 
ders of magnitude more than the second one as deter- 
mined by the transformation matrix A4 ti> j V . At a bias of 
about 0.5V, all vibrational modes of the molecular junc- 
tion are excited. The mode that contributes the most 
to the heating corresponds to an in-plane "breathing" of 
the central carbon ring (not shown). |To| 

So far we have assumed the junction is thermally iso- 
lated. However, the majority of the heat generated in 
the junction is actually transferred to the electrodes. We 
estimate the thermal conductance following the approach 
of Patton and Geller. [ll| We assume the junction to be 
a weak thermal link with a given stiffness K, which we 
evaluate from first principles. ^3] We then estimate the 
thermal current into the electrodes via elastic phonon 
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FIG. 4: Local temperature as a function of bias for a benzene- 
dithiolate molecular junction (dashed line) and a gold point 
contact (solid line) due to the equilibrium between local heat- 
ing and heat dissipation into the electrodes at T e = K. 



scattering as [l"]| 
47t.fr 2 

hh — r — 



deeN L (e) N R (e) [n t (e) - n R (e)] , (13) 



where is the Bose-Einstein distribution function, 

Nj^fm (e) is the spectral density of states at the left 
(right) electrode surface, and K = Trd 2 Y/(Al). Y is the 
Young's modulus of the junction, and d (I) is its diam- 
eter (length). [llj We note that the temperature pro- 
file along a nanojunction is nearly constant and almost 
equal to the average temperature of the thermal reser- 
voirs. |l3l FLU Hil Hq The thermal current from the junc- 
tion with temperature T w dissipated to the left electrode 
with temperature T e is thus equivalent to the thermal 
current of a weak thermal link between reservoirs with 
temperatures T e and 2T W —T e . Analogously for the ther- 
mal current into the right electrode. The effective lo- 
cal temperature, resulting from the equilibrium between 
local heating and heat dissipation into the electrodes, 
is plotted in Fig. 01 at T e = K. Values of the local 
temperature for the gold contact are in good agreement 
with experimental results 6] and previous theoretical es- 
timates, pj It is evident from Fig. 0] that the local tem- 
perature is considerably reduced by heat transfer into the 
electrodes even for relatively large biases (cf . Fig. [5J ■ It 
is also evident that at any given bias the single-molecule 
junction heats up less than the gold point contact. This 
is due to the larger stiffness of the molecule and its larger 



resistance to electrical current. 
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